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Abstract 

We investigate 0(4) textures in a background with a positive cosmological 
constant. We find static solutions which co-move with the expanding back- 
ground. There exists a solution in which the scalar field is regular at the 
horizon. This solution has a noninteger winding number smaller than one. 
There also exist solutions in which scalar-field derivatives are singular at the 
horizon. Such solutions can complete one winding within the horizon. If the 
winding number is larger than some critical value, static solutions including 
the regular one are unstable under perturbations. 
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I. INTRODUCTION 



Phase transitions associated with the symmetry breaking in the early universe could 
produce topological defects. Those are core defects such as domain walls, cosmic strings, 
and monopoles, or textures |1J]. 

Global textures are produced when a continuous global symmetry G is broken to a group 
H. The resulting homotopy group ttd(G/H) is nontrivial, and the vacuum manifold is S D ; 
for G=0(N) models, if =0(N-1) and D — N — 1. The mapping from the physical space to 
the vacuum manifold is R D — > S D . 

The Lagrangian describing such texture models is 

£ = --d^ a d fl $ a --($ a $ a -T] 2 ) 2 , (1) 
2 4 

where a = 1,...,N. This Lagrangian has O(N) symmetry. After the phase transition, the 
scalar field acquires the value in the vacuum manifold, $ a $ a = t] 2 , of which symmetry is 
O(N-l). 

The energy of this model in D spatial dimensions is given by 
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E = I d D x 



2 (V$ a ) 2 + 1/($ a ) 



h + h ■ (2) 



By transforming $ a (x) under the rescaling x — > ax, the energy configuration becomes 

E a = a 2 ~ D h + a- D I 2 . (3) 

For a global texture, the potential part of the energy vanishes, I2 = 0, since the scalar field 
sits on the vacuum manifold in the entire space R D . Equation (^) tells us that the energy 
is a function of the rescaling factor a. The configuration relaxes to acquire the minimum 
energy. Therefore, 0(2) texture (D = 1) expands (a decreasing), and 0(4) texture (D = 3) 
collapses since there is no stable minima in E a with respect to (w.r.t.) a. Only 0(3) texture 
(D = 2) is stable. This is the well-known Derrick's theorem 0. 

Collapsing 0(4) textures produced during the Grand Unified phase transitions attracted 
much interest in the early 1990's due to their potential cosmological significance ||-|6|]. As 



the universe expands and cools down, the texture knots shrink and radiate away into Gold- 
stone bosons. The texture collapse occurs at the speed of light, and correlates the scalar 
field on larger and larger scales. Quickly the correlation length grows to the horizon scale, 
and expands with it. As the texture collapses, the gradient energy of the scalar field gets 
concentrated at the center of the texture. This energy becomes large enough for the scalar 
field to overcome the potential barrier, and to leave the vacuum manifold. The texture un- 
winds in the end. During the collapsing texture leaves a signature in the structure formation 
of the universe. 

Textures do not require a complete winding for collapse. Above some critical noninteger 
winding number ~ 0.5, textures collapse |7|,§]. The stability of textures in an expanding 
background has also been investigated for the radiation-dominated, matter-dominated and 
inflation era IP HTOfl - 

A non-vanishing cosmological constant has attracted increasing attention as a source for 
an accelerating universe. Supernova observation data show that the universe is accelerating 
rather than decelerating. The simplest source for such expansion is a cosmological constant. 

If it existed, the story of texture stability could be altered. With a cosmological constant 
the universe undergoes de Sitter expansion. 0(4) textures formed in such a background could 
slow down their collapse, or could even stabilize themselves to co-move with the background. 

In this paper, we investigate 0(4) textures with a cosmological constant. We study both 
the scalar and the gravitational field configurations. We solve the field equations numerically, 
and seek the static solutions. We also check the stability of them. 

When there is a positive cosmological constant, a horizon develops in the geometry. There 
exists a static texture solution which is regular across the horizon. However, the scalar field 
of this solution does not complete its winding beyond the horizon; the winding number is 
nonintegral ~ 0.88. There also exist static solutions in which scalar-field derivatives are 
singular at the horizon; however, the scalar field itself is "finite" at the horizon. The value 
of the scalar field at the horizon can be either higher than that of the regular solution, 
or lower. For this type of solution, the texture can even complete one winding within the 
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horizon. Both of the regular and the singular solutions are physical. 

The regular solution is unstable under a small perturbation. The singular solutions of 
the high scalar field at the horizon are also unstable. The others look stable. 

In Sec. II, we present the model and field equations. In Sec. Ill, we study texture solutions 
in pure de Sitter space. In Sec. IV, we include the texture self-gravity. In Sec. V, we check 
the stability of the solutions. In Sec. VI, we discuss the texture evolution. We conclude in 
Sec. VII. 



II. THE MODEL 

The action describing an 0(4) texture with a cosmological constant is 

'TZ-2K 1 A 



-d u <5> a dv$ a - - n 2 ' 2 



16nG 2 " 4 



(4) 



where $ a is the scalar field with a = 1,...,4, rj is the symmetry breaking scale, and A is 
the cosmological constant. We are mainly interested in the case of A > 0. We assume the 
spherical symmetry and the general static form of the metric is 

ds 2 = -B(r)dt 2 + A(r)dr 2 + r 2 dn 2 2 . (5) 

We consider a scalar field everywhere in space to reside in the vacuum manifold, $ a $ a = 
rf. The scalar field is of the form, 

$ a = i](sm x sin 9 cos (f>, sin \ sin 9 sin 0, sin \ cos 9, cos x) , (6) 

where 9 and are the usual spherical coordinates, x is a function of the radial coordinate, 
x(r). In this case, the nonlinear cr-model approximation is applied, and A in the potential 
is treated as a Lagrangian multiplier. The scalar field equation, then, is 

The energy-momentum tensor is given by 
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t, u = d^ a d u $ a - l - g ^d a § a d a <$> a . (8) 

With the given metric fl5|), the Einstein's equation leads to 

-«i4f^^^ ! (^^) + A, (9) 



A \ At r 2 J ' \2A r 2 

B' A-l\ o ^ 2 / x' 2 sinV 
A\ B~r + — =87rGV \'2A + — 2 



~ 9i ~ A 



B" 1 (B'V A'B' B' A 



2B 4\B J AAB 2Br 2Ar 
and the scalar field equation leads to 



SttG^+A, (11) 



y" 1 ( B' A' 2\ . sin2 X n 

In the field equations, A can be absorbed into the radial coordinate by rescaling. Then, 
the equations are free of A, and the solutions are as well. This implies that the solutions 
look self-similar one another for different A scales. The only physical difference is that the 
configuration stretches along the radial direction by the A scale. 

III. SOLUTIONS IN PURE DE SITTER SPACE 

In this section, we investigate texture configurations in pure de Sitter space without 
including the self-gravity of the texture field. In this case, solutions are more tractable 
analytically as well as numerically. The background geometry is 

ds 2 = -(1 - -r 2 )dt 2 + 1 dr 2 + r 2 dQ 2 2 . (13) 

3 1 — ^r 2 

The scalar field equation with this metric is 

20-2^) = 

X R(l- R 2 ) X R 2 (l -R 2 ) K } 

Here, we rescaled the radial coordinate by R = (A/3) 1 ^ 2 r, and the prime denotes differenti- 
ation w.r.t. R. The horizon is located at R = 1. This equation is solved only numerically. 
The numerical calculation is impossible at the horizon, R — 1, where the equation becomes 



singular. However, if we fix the value of the scalar field at the horizon as a boundary con- 
dition, we can avoid the difficulty. The equation is numerically solved by the relaxation 
method after fixing two-point boundary values at R = and R — 1. We locate the center 
of texture at the origin, x(0) = 0. 

First, let us discuss the regular solution. If x is finite up to the second derivative at the 
horizon, from Eq. (|14|) we have a regularity condition for x a t the horizon as 

2x' h + sm2 Xh = 0, (15) 

where the subscript h denotes the value at the horizon. We found numerically that this 
condition is satisfied when Xh — 0.641927T. If we proceed numerical calculation outside the 
horizon with this value, the scalar field reaches x{R ~^ °°) ~ 0.720817T asymptotically. The 
scalar field of this regular configuration does not complete winding in the vacuum manifold 
even at infinity. The solution is plotted in Fig.|I|(a). 

Besides the regular solution, there exist solutions in which x" & t the horizon is infinite. 
From Eq. (0), when x" is infinite, unfortunately there would not be a definite information 
of x an d x' & t the horizon; x' can be either finite, or infinite at the horizon. However, we 
observe that x' approaches a very large value close to the horizon. This indicates that x' 
may be infinite there — we checked that at least x" (R ~ l) s type of solution does not 
exist, where s < 0. 

Such solutions exist for Xh ~ 0.4457T and 0.657T < Xh ~ 1.17T. However, close to Xh ~ 
and ~ 0.657T the derivative x'h * s n °t very large though. Other values of Xh than those 
for regular and singular solutions, field configurations are not so conventional; they do not 
monotonically increase, but decrease close to the horizon. The singular solutions are plotted 
in Fig.Q(b). 

Note that in flat spacetime as well as in self-gravitating spacetime without a cosmological 
constant, there do not exist static configurations above some critical noninteger winding 
number smaller than one. However, with a cosmological constant, we have static solutions 
up to Xh ~ l.lvr. For Xh = tt, the texture completes one winding at the horizon, and for 
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7r < Xh ~ 1.17T, it does even inside the horizon. We will see that this is also true for the case 
including the self-gravity of the texture field in the next section. This result agrees with the 
expectation we mentioned in the Introduction. Collapse is balanced by the expansion of the 
background. 

For a given singular configuration inside the horizon, the matching solution outside must 



exist. However, as the field equation ([L4]) does not provide any information about the 
functional form, or specific values of the field derivatives at the horizon, we cannot continue 
numerical calculations to the region outside of the horizon. 

Considering solutions only inside the horizon is justified as following. When the texture 
is formed, the scalar field acquires values in the vacuum manifold. Different regions in space 
take different values of the field. The correlation length scale within which the region takes 
the same value, cannot exceed the cosmological horizon scale at the moment of formation. 
Usually the formation of one defect is completed within the cosmological horizon. Other 
defects are formed in the neighborhood. 

In the present case, if we include the contribution of other matters (for example, the 
texture field itself) to the energy density of the universe, the cosmological horizon deter- 
mined by the total energy density is always smaller than the de Sitter horizon. Therefore, 
considering solutions only inside the de Sitter horizon is perfectly guaranteed. 

Although singular solutions would be accepted by the above justification, we need to 
comment some trouble associated with them. First, if x' is infinite at the horizon, the 
curvature diverges. If we include the self-gravity of the texture, the horizon is not only a 
coordinate singularity, but becomes a curvature singularity. If there is a curvature singular- 
ity, the signal emitted from it can propagate to the inner region. Observers off the center 
of the texture eventually will see the signal. This trouble may be resolved if we put the 
causal boundary to the cosmological horizon which is inside the de Sitter horizon. Hence, 
the texture configuration is cut by neighboring textures before it develops a singularity. 

As we will see in Sec. V, the static texture configurations are unstable under small 
perturbations. The scalar field may evolve in such a way as to remove the singular behavior 



at the horizon. 

In addition, when the field gradient is very large locally, the scalar field acquires enough 
energy to overcome the potential barrier. Therefore, the nonlinear a-model approximation 
is not valid any longer. 

Even though the solutions obtained here are static inside the horizon, they co-move with 
the expanding background. The configuration will be self-similar as time passes unless there 
are perturbations. 



In this section, we include the effect of the texture self-gravity. For textures of GUT 
symmetry-breaking scale, 77 ~ 10~ 3 m p ~ 10 16 GeV, the effect is very small. Therefore, the 
results are very similar to those in the previous section. 

The location of the horizon is displaced a little bit from that of pure de Sitter space. 
Numerically we can perform calculations close to the horizon where 1/A « 10~^ 5 \ Be- 
yond that point, the values in calculation exceed the numerical limit. Therefore, we obtain 
numerical solutions up to that point. 

Let us investigate the regular solution at the horizon. If we define Q = 1/A, the scalar 
field equation ( |12D becomes 



All physical quantities in the above equation are regular at the horizon. We perform Taylor 
expansion for each quantity about the horizon, R = R h + x, where B(R h ) = Q(Rh) = 0, 



IV. SELF-GRAVITATING SOLUTIONS 




(16) 



B 



b\x + b 2 x 2 + b 3 x 3 + ... , 



(17) 



Q 



q x x + q 2 x 2 + q 3 x 3 + ... , 



(18) 



X 



fo + fiX + f 2 x 2 + f 3 x 3 + .... 



(19) 



We plug these quantities in the field equation flTB|), and expand it in the power of x. Equating 
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the coefficient of each term to be zero gives relations among coefficients in the expansion 
series above. The first term, 1/x, in the expansion gives 

In other words, 

Xh+^L 2 Xh = 0. (21) 

Note that this condition reduces to the regularity condition ( |i"5"D in de Sitter space for Rh = 1 
and A = 1/(1 —R 2 ). We numerically search the solution which satisfies the above condition 
close to the horizon. We found that the solution exists when x{Rb) — 0.642950227T for 
rj = 0.02m p ,[| where Rb is the location of the last mesh of our numerical calculation. This 
is very close to Xh for the pure de Sitter solution. The solution is plotted in Fig.|2|. For this 
solution Rb ~ 0.99535 and l/A(Rb) ~ 8.7 x 10 -5 . Although we are not able to perform 
numerical calculations beyond the horizon, we believe that the configuration is also very 
similar to the pure de Sitter one there. 

We have also singular solutions similar to the de Sitter ones. For the singular ones, we 
performed the calculation up to Rb ~ 0.99483 where \/A{R^) ~ 1.0 x 10~ 4 . For the texture 
with one winding number, we imposed the boundary value of the scalar field, = 
0.92827T. Since x increases very rapidly close to the horizon, it is reasonable to impose this 
boundary value at R = Rb estimated from the de Sitter solution for Xh = This solution is 
plotted in FigfJ as well as the other singular solution for x{Rb) — 0.531497T. Similar to the 
de Sitter solutions, there also exist unconventional solutions for 0.531497T < x(-Rfc) ~ 0.647T. 
This range got narrowed compared with the de Sitter one. 

The gravitational fields of the singular solution for x{Rb) = — 0.92827T are plotted in 
FigS They are not much different from the pure de Sitter solutions. 



Throughout the paper, we show the numerical results for rj = 0.02m p in order to see the gravi- 
tational effect clearly. 
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As it was discussed in the previous section, the singular behavior of the derivatives of the 
scalar field at the horizon turns the horizon to be singular. From the regularity condition at 
the horizon, x" is infinite. We cannot determine the definite value of %' there. It could be 
either finite, or infinite. From our numerical calculations, however, it is obvious that \' is 
a very large number. This indicates that the curvature at the horizon diverges. The Ricci 
scalar evaluated from the energy-momentum tensor depends on x', 

n = 16nGr]2 (iz + ^) + 4A • (22) 

We plot the Ricci scalar for the singular solution in Fig.[|. It diverges steeply approaching 
the horizon. The other scalar curvatures evaluated from the metric contains x' ■ However, 
it is not so clear if they involve the terms of x" which is clearly divergent. 

The horizon in pure de Sitter space is an observer-dependent property. To a given 
observer A, it is A's future and past event horizon beyond which events are causally dis- 
connected with A. Observers at different locations in space have their own horizon which is 
different from one another. If there is a texture centered at A's location, the spacetime is no 
longer homogeneous. In this situation, observers have their own horizon which is different 
from the one in pure de Sitter space. Their horizon is still a local property in terms of 
causality. However, the curvature singularity located at A's horizon is a global property. 
Every observer should recognize this singularity if it is located inside his/her own horizon. 
What makes the horizon, a priori, a local property of A's be so special (singular) is that 
the texture is "centered" at A's location. The way to circumvent the trouble caused by the 
singularity was discussed in the previous section. 

It is not unusual that a naked singularity exists in the geometries of topological defects. 
For example, global cosmic strings and supermassive gauge cosmic strings possess a singu- 



larity in the spacetime [TT|,[T2|. The difference is that the singularity in the string geometry 
is accompanied with the regular scalar field configuration. In the case of texture, the sin- 
gularity arises only when the scalar field behaves singularly at the horizon; the horizon is 
regular if the scalar field is regular there. 
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V. STABILITY 



In previous sections, we found out static texture solutions co-moving with the expanding 



background. In Refs. |8|-10|, it was investigated that textures are very unstable even in 
expanding backgrounds. Textures collapse, or expand further. It is worthwhile to check if 
our static solutions are stable, or not. We adopt the linear-perturbation method used, for 



example, in Ref. [I3"| where authors check the stability of black holes. This method finds 



if there is an unstable mode in perturbations, but does not describe the evolution precisely 
after the perturbations. 

Let us expand the gravitational fields and the scalar field as 

B(t, r) = B (r) + eB^t, r) + 0(e 2 ) , (23) 
A(t,r) = A (r) + eA^r) + 0(e 2 ) , (24) 
X (t, r) = xo(r) + exi(t, r) + 0{e 2 ) , (25) 

where B (r), A (r), Xo( r ) are the static solutions we obtained in the previous section, and 
e is a small number. We expand the time-dependent field equations in the power of e, and 
keep the terms only up to the first-order in e. The time-dependent Einstein equations are 



A \Ar r 2 1 \2B 2A r 2 ' 



_ G ;; = U-*+A=l) = ftrfV ( _ *1 _ ^ + ) + A , (27) 
r A\ Br r 2 1 \ 2B 2A r 2 ' 



G tr = ^=8nGri 2 xx'i (28) 



and the scalar field equation is 

B B\2B 2A X A A\2B 2A + r X r 2 



We plug expansions (|23)-(|2~5) in the above field equations. The first-order terms in e gives 
equations as following. 
From the G tr equation, 
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Al 



A r 



i • 



(30) 



of which solution is 



Air 



(31) 



where f(r) is an integration constant w.r.t. time. This solution satisfies G\ equation and we 
can fix the integration constant, fir) = 0, as a boundary condition. This gives the unique 
solution to Eq. fl3lD . This boundary condition is appropriate because it is consistent with 
reasonable initial conditions, A(t = 0,r) = A (r) and x(t = 0,r) = Xo( r )- This means that 
perturbations are zero initially, Ai(0,r) = and Xi(0, r) = 0, but their velocities are not. 
The G\ + G r r equation gives 



1 (B{ A[ B'^B, A! A 1 
A I S A Bl Al 



2 sill2 X0 

lo7rG?7 Xi 



(32) 



The scalar field equation gives 
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+ ; 



2 cos 2%o 



A 2 , 



, . ^isin2xo n 

*o + 1 — = 

Aq r l 



(33) 



Using Eqs. ([H|) and (|32]) , the perturbation terms B\, A 1 , B[, and A[ are nicely eliminated 



in the scalar field equation (|3~3|), 



4> 



■B-Xi - Xi + a(r)x[ + &(r)xi = , 

DO 



(34) 



where 



a(r) 



B' 



+ 



A' 



b(r) = Ao 



2B 2A 
2 cos 2xo 



2 

r 



+ 87rG V 2 x' r 



2 sin 2xo 



1 



B' 



An 



~ Xo 



(35) 
(36) 



r 2 A V2B 2 A, 

Therefore, the perturbation equation for xi involves only the quantities obtained from static 
solutions. Let us take the perturbed scalar field in the form, 
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Xi{t,r) = (p(r)e 



iujt 



(37) 



then the equation for the amplitude (p(r) becomes 



B 



^-[—f"{r) + a(r)(f'(r) + b{r)<p{r)\ = u 2 ip(r) . 
Aq 



(38) 



This is an eigenvalue equation with the eigenvalue u 2 . We rescale the radial coordinate and 
amplitude function by 



An 



ij){z) = N 



o V 25 

(p(r)r 



dr , 



(39) 
(40) 



where N is a normalization constant. This rescaling puts the equation into the non 
relativistic Schrodinger-type equation in spherical coordinates 

Id 2 Id 



2 dz 2 z dz 



+ V(z) 



ip(z) = 10 ip(z) 



(41) 



where the potential is 



V[z(r)] = ^ 



" 1 (B' Q A'A . 



(42) 



The advantage of deriving the equation into this form is that the shape of the potential 
V provides insight for the eigenvalue spectrum. If V is a potential well, there exists a 
discrete eigenvalue spectrum. If the minimum of V is negative, there possibly exists negative 
eigenvalues, so u can be imaginary. In such a case, the perturbed field has an exponentially 
increasing, or decreasing solution depending on the boundary condition we impose. Even if 
there exists only one negative eigenvalue of u 2 , the system becomes completely unstable. 

Now with the given equation above, let us examine the stability of the solutions obtained 
in the previous section. First, for the regular solution, the potential V is plotted in Fig.|5]. It is 
a potential well with a negative minimum. We numerically solve the eigenvalue equation (p5[) 
rather than Eq. fljip. As the potential V diverges at origin, we set the boundary condition, 
■0(0) = 0, so ip(0) = 0. For the regular solution at the horizon, x'o( r h) is finite, and B (rh) = 
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l/Afrh) = 0. Therefore, it is easy to see <p(rh) = 0. With these two boundary conditions for 
<p(r), we found a negative eigenvalue, u 2 = — 0.7033105A. The corresponding eigenfunction 
ip(r) is plotted in Fig.|5|. Our regular solution is unstable under small perturbations. 

Now, let us examine singular solutions. We are particularly interested in the solution 
with Xo{ r h) ~ 71 , which completes one winding in the region under consideration. The 
potential V is plotted in Fig|| The potential diverges at the origin similar to that of the 
regular one. And, as we approach the horizon, r — > (z — > oo),0the potential also diverges. 
Therefore, we impose boundary conditions, ^(0) = ip( z ~~ > °°) — 0? which, in turn, gives 
y?(0) = (p{rh) = 0. Numerical calculation shows that there exists a negative eigenvalue, 
uj 2 = — 0.207815A. The corresponding eigenfunction is plotted in Fig.|6|. The static singular 
solution with Xo( r h) ~ vr is unstable. 

For the other singular solutions with < Xo( r b) ~ 0.531497T, the potential V has a 
negative minimum also. However, it is very close to zero. There does not necessarily exist a 
negative eigenvalue in such a case. In particular, if the symmetry-breaking scale rj is raised, 
the minimum value of V is raised above zero. Therefore, we conclude that this type of 
singular solutions is stable. 

Even though we are not so convinced to study the unconventional solutions with 
0.531497T < Xo( r b) ~ 0.647T, the potential V for those solutions has a deep negative mini- 
mum, and the solutions are unstable. Therefore, in the rest of this paper, we will consider 
that textures are stable only for Xb ~ 0.531497T. 

As a result, the static regular solution and singular solutions with a high winding number 



2 Even though finite z region is shown in the plot, z goes to infinity as we approach the horizon. 
This is evident if we consider the pure de Sitter case where analytic form of radial rescale is 
possibly obtained, z = (r/j/\/2) tanh~ 1 (^/r/r/ l ). We insist that the situation be similar in the 
self-gravitating case. 
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up to one (or, a little bit higher than one) are unstable under perturbations.^ Depending 
on the boundary conditions imposed at the moment of perturbations, textures collapse, or 
expand w.r.t. the expanding background. 

The perturbed expansion proceeds behind the cosmological horizon, so it is not very 
interesting cosmologically. The collapse will be discussed in the next section. 

VI. COSMOLOGICAL EVOLUTION 

In this section, we qualitatively discuss the evolution of the configurations we obtained. 
As we mentioned earlier, the configurations are self-similar regardless of the A scale. Impos- 
ing a definite A value fixes the physical length scale. 

If A is very large and dominates the energy density of the universe, the situation corre- 
sponds to inflation. The large A plays the role of the vacuum energy given by the inflaton 
field. In this case, the cosmological horizon scale determined by the total energy of the 
universe is comparable to the horizon scale ~ (A/3) -1 / 2 of our solution. Then the whole 
configuration we obtained up to the horizon is applicable. 

Unstable textures collapse during inflation. Even if the collapse reaches a relativistic 
speed, however, it would not be able to compensate the inflationary expansion. Textures 
survive for the inflation era. After inflation ends, textures will collapse to the origin very 
quickly, and unwind themselves in the end. 

For the rest of the history of the universe except present, A remains very small compared 
with the total energy. The cosmological horizon scale is much smaller than r^. If textures are 
formed at this stage, the causality bound prohibits the formation of a large texture extended 
to r/j. Neighboring textures are formed before that. Even though such a large texture had 



3 The winding number is evaluated as W = Xb/ft — (sin 2xb)/^, where Xb is the value of the 
scalar field at the boundary under consideration. In this work, however, we shall mostly use Xb f° r 
numerical values rather than elaborating to calculate the corresponding W. 
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been formed with a high winding number up to one, the observers would see only a part of 
it, so an incomplete winding in the universe. After textures start to collapse, they will see 
the scalar field with higher values crossing into the boundary of the universe. In the end, 
the observers who initially saw a part of the winding, will see the complete winding of the 
field. 

We found that the static textures for Xb ~ 0.531497T are unstable. Even though the 
universe with a cosmo logical constant is topologically different from those without it, (for 
example, there is no static solution for high winding numbers in the others) collapsing 
behavior is very similar. In Ref. ||, it was investigated that textures collapse for winding 
number W > W c ~ 0.6 in the flat, radiation, and matter background. Depending on the 
initial configuration imposed, the collapsing time varies. It can be explained as following. 
For W > W c , there does not exist a static solution in such backgrounds. Imposing the initial 
configuration already implies a certain amount of perturbation. The higher W is, the larger 
are the perturbations. For our case, the static configuration will start to collapse slowly. 
However, the collapse will reach a relativistic speed in the end. 

Textures with \b ~ 0.531497T are stable according to our investigation. As it was pointed 
out in Ref. ||, they behave as a local fluctuation, and will decay into Goldstone bosons. 

VII. CONCLUSIONS 

In this work, we investigated 0(4) textures with a positive cosmo logical constant. It has 
been previously investigated that 0(4) textures collapse in the flat, self-gravitating, radia- 
tion, and matter background. There was no static solution found out in those backgrounds, 
especially for high winding number textures. However, with a cosmological constant, we 
found out that 0(4) textures have static configurations for all winding numbers up to slightly 
higher than one. De Sitter expansion of the background prevents the texture from collaps- 
ing. However, the static configuration is unstable under small perturbations if the winding 
number is larger than some critical value, Xb ~ 0.53 1497T (for rj = 0.02m p ). 
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The textures with such high winding numbers will collapse after the perturbations are 
applied. Initially the collapse will occur slowly, but it will reach a relativistic speed in the 
end. At the final stage, the field gradient energy will be concentrated about the center, and 
textures will unwind. Meanwhile, textures will leave a signature in the large-scale structure 
formation of the universe. The resulting signature will be similar to that of the textures in 
the other backgrounds. 

For Xh ~ 0.53 1497T, textures are stable. They will decay into Goldstone bosons in the 
end. 

The static configurations we obtained fit the region inside the horizon. They co-move 
with the expanding background while they are stable. The scalar field is definite everywhere 
within the horizon. For most of the configurations except the regular one, the static solutions 
have singular field derivatives at the horizon. This singular behavior turns the horizon to 
become a curvature singularity. However, this is not a matter so much because the horizon 
is always located outside the cosmological causal bound. 

The detailed evolution based on solving time-dependent field equations was not inves- 
tigated in this work. It will be interesting to perform numerical simulations for such a 
evolution. In this case, we can extend the model to a quintessence field background of which 
equation of state is P = wp (—1 < w < 0). For such a field, there would not exist a static 
description of the spacetime. Therefore, numerical simulation is very useful to investigate 
the evolution. The quintessence field may affect the texture evolution more than the cos- 
mological constant, because its energy density could be relatively high at the early stage of 
the universe. 
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FIGURES 




FIG. 1. Plot of the scalar field in pure de Sitter space, (a) Regular solution. At the horizon 
the scalar field has Xh — 0.641927T. It approaches asymptotically \ ~ 0.720817T. (b) Singular 
solutions for \h = tt and Xh = 0.4457T. We also show an unconventional solution for Xh = 0.57T in 
which the scalar field decreases close to the horizon. 
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FIG. 2. Plot of the self-gravitating scalar field for r] = 0.02m p . The middle one is the regular 
solution with x(i?&) — 0.642950227T. The others are singular solutions for x(Rb) = 0.9282-7T and 
x(R b ) = 0.53149tt. 
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FIG. 3. Plot of the gravitational field for 77 = 0.02m p and \{Rb) = 0.9282tt. B and l/A are 
very close to each other as well as to the pure de Sitter field 1 — R 2 . 
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FIG. 4. Plot of the Ricci scalar for the singular solution with 7] = 0.02m p and x{Rb) = 0. 
It looks divergent as it approaches the horizon. The dotted line is for the regular solution. 
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FIG. 5. Plot of the potential V for the perturbation field of the regular solution with 
rj = 0.02m p . There exists an unstable mode with the eigenvalue to 2 = — 0.7033105A. The cor- 
responding eigenfunction tp is also plotted. The normalization of tp was not taken into account. 
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FIG. 6. Plot of the potential V for the perturbation field of the singular solution with 
r] = 0.02m p and x(Rb) = 0.9282-7T. There exists an unstable mode with the eigenvalue 
lj 2 = — 0.207815A. The corresponding eigenfunction ijj is also plotted. The normalization of tp 
was not taken into account. 
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